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We present an unquenched lattice calculation for the — B" transition amplitude. The cal- 
culation, carried out at an inverse lattice spacing 1/a — 2.22(4) GeV, incorporates two flavors of 
dynamical quarks described by the fully 0(a)-improved Wilson fermion action and heavy quarks 
described by NRQCD. A particular attention is paid to the uncertainty that arises from the chiral 
extrapolation, especially the effect of pion loops, for light quarks, which we find could be sizable 
for the leptonic decay of the Bd meson, whereas it is small for the Bg meson and the B param- 
eters. We obtain fs^ = 191(10)(t^^) MeV, /s,//b^ = 1.13(3) (t'l), BsAmb) = 0.836(27)(1^|), 
BbJBb^ = 1.017(16)(+i7), C = 1.14(3)(1"), where the first error is statistical, and the second is 
systematic, including uncertainties due to chiral extrapolation, finite lattice spacing, heavy quark 
expansion and perturbative operator matching. 

PACS numbers: 12.38.Gc, 12.39.Fo, 12.39.Hg, 13.20.He, 14.40.Nd 



The unitarity test for the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix entered a new era with the 
BaBar and Belle measurements of the angle (pi P, 0- 
The test requires the determination of the other an- 
gles and the sides of the unitarity triangle, the precision 
of the latter being limited by uncertainties in hadronic 
matrix elements. Lattice QCD in principle offers a 
model-independent calculation of such matrix elements. 
Those matrix elements provided by lattice calculations 
to date, however, are based on the quenched approxima- 
tion, blindly hoping that quenching does not introduce 
large errors. 

Simulations including creation and annihilation of a 
quark anti-quark pair in the vacuum have become feasible 
only recently. In this letter we present an unquenched lat- 
tice calculation of the hadronic matrix elements appear- 
ing in the — mixing amplitude, which is needed 
in the determination of the CKM matrix element \Vtd\ 
from the mass difference AM^ j^. The matrix element 
is parametrized as (i?°|67^(l - 75)9^7a'(1 - 75)g|Sg) = 
f /s BB^Mg , where /b, is the B meson decay constant 
and Bb^ denotes the B parameter (q represents d or s 
quark). Our prime interest is to calculate Bb^, as un- 
quenched calculations of Jb are already available 0. 
We include, however, the calculation of Jb to provide a 
consistent set of the matrix element for B'^ — i?" mixing. 

With current lattice calculations, systematic uncer- 
tainties due to the discretization error and the pertur- 
bative matching between continuum and lattice oper- 
ators amount to 10-20%. One may improve the ac- 



curacy of \Vtd\ by studying the ratio AMg/ AMd if 
Bg — Bs mixing is measured. The relevant quantity is 
£. = (/b, V ^B, )/{fBa \jBBa ) , where many of the sys- 
tematic uncertainties cancel in the lattice calculation. 

A remaining major uncertainty arises from the chiral 
extrapolation of the lattice simulation which is made with 
relatively heavy dynamical quarks. One may resort to 
chiral perturbation theory (ChPT) as a theoretical guide 
for the extrapolation. The problem is that the currently 
available lattice data do not show the logarithmic be- 
havior expected from long-distance pion loops in ChPT 
[3 . It is in the scope of the present work to discuss the 
uncertainty in the the matrix elements associated from 
the chiral extrapolation in the absence of the observable 
logarithmic behavior. 

The calculation is carried out on the unquenched gauge 
configurations generated at /? = 5.2 on a 20'' x 48 lattice. 
Two flavors of dynamical quarks for the u and d quarks 
are simulated at five values of quark mass in the range 
(0.7 — 2.9)ms with the physical strange quark mass. 
This corresponds to the pseudoscalar to vector mass ra- 
tio of 0.6-0.8. The hopping parameter chosen is i^^ea 
= 0.1340, 0.1343, 0.1346, 0.1350, and 0.1355. For each 
sea quark mass, 1,200 configurations are accumulated for 
measurements from 12,000 HMC trajectories separated 
by 10 trajectories. The lattice spacing we adopt is de- 
termined from p meson mass and equals 2.Tl(A) GeV 
after the extrapolation to the chiral limit. This value is 
consistent with 2.19(15) GeV from the Sommer scale ro 
assuming the physical value of 0.49 fm, and 2.25(5) GeV 
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FIG. 1: Chiral extrapolation of (filled circles) and 

(open squares) . The quadratic extrapolation is shown by solid 
lines, while the fits with the hard cutoff chiral logarithm are 
shown for \i = 300 (dotted curve), 500 (thin dashed curve) 
and cxD (thick dashed curve) MeV. Quenched results are also 
shown (triangles). 



from Jk (with an additional 0(5%) error from the per- 
turbative matching). This suggests that the large width 
of p may not seriously affect the chiral extrapolation of 
the p meson mass. Other details of the simulation are de- 
scribed in [^. We adopt the lattice NRQCD formahsm 
0, 111 for heavy quarks. The non-perturbatively 0(a)- 
improved Wilson action is employed for both valence 
and sea light quarks. 

We take five values of heavy quark mass mq {artiQ = 
1.3, 2.1, 3.0, 5.0 and 10.0) to cover mQ= 3-20 GeV for the 
NRQCD action that contains all corrections of the order 
l/niQ. We take the valence light quark mass set equal to 
the sea quark mass and then extrapolate to the physical 
u- and d-quark masses, unlike the 'partially quenched 
analysis' often adopted in the literature. The strange 
quark is treated in the quenched approximation. There is 
an uncertainty in the determination of the strange quark 
mass depending upon which strange hadron is used as 
input. We take three values for the hopping parameter 
K, ^ 0.13465, 0.13468 and 0.13491 and interpolate to 
Ks{K) — 0.13486(3) (for the K meson as input) and 
K,{(t)) = 0.13471(4) (for the meson as input) The 
method to calculate Jb and the B parameter follows our 
previous studies in the quenched approximation |lOlllll| . 

Figure ^ shows the chiral extrapolation of the decay 
constants /b, expressed in terms of $/^^ = fs^ \JM 
as a function of the pion mass squared. For Jb,, the 
result is shown at Ks=Q.134:65. In order to absorb the 
change of effective lattice scale that varies with Kgea at a 
fixed bare coupling /3, both axes are normalized with rg 
determined from the heavy quark potential at each sea 
quark mass. The heavy quark mass is interpolated to the 
b quark using the lattice data. 

Open triangles show quenched results obtained at a 



similar lattice spacing 1/a = 1.83(2) GeV {f3 = 6.0) with 
the non-perturbatively improved Wilson quark action. 
Our observation that they lie close to the unquenched 
data (filled circles) implies that the B meson decay con- 
stant takes a similar value in quenched and two-flavor 
QCD if the scale is normalized by tq. With the more 
conventional normalization of using the p mass, however, 
the unquenched values are higher by about 20% (see, e.g., 
4]), as is seen by the fact that rouip = 1.91(2) for Nf 
— 2 while it is 2.20(3) on the quenched lattice This 
is understood as systematic errors of the quenched ap- 
proximation, with which the determination of the lattice 
scale depends on which physical quantity is the input. 
With the dynamical quarks these errors are significantly 
reduced, leading to a convergent determination of the 
lattice scale. 

The solid line represents a linear plus quadratic fit in 
(roTTiTr)^, which describes the lattice data well. This fit, 
however, does not contain the chiral logarithmic term 
which is predicted by ChPT (IJJ: 
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where terms regular in are omitted, and the coupling 
q is the B*B'K interaction in ChPT. Lattice calculations 
|l3l Il3| give a value consistent with the empirical one 
measured for D* Dn decay, g = 0.59 ±0.01 ±0.07 
Although there is an uncertainty in translating the value 
at the I?'^*-' meson to that at the infinitely heavy quark 
mass (where the heavy-light ChPT is formulated), we 
take g = 0.6 to estimate the effect of the chiral logarithm. 

Let us here consider a simpler case. For the pion de- 
cay constant the one-loop logarithmic term is controlled 
by the number of dynamically active quark fiavors Nf as 

— 2 (4^/)2 In -jjT ; no uncertam parameters such as g are 
involved. Thus, the test for the presence of the chiral 
logarithm is less ambiguous j^. Our high statistics un- 
quenched data, shown in Figure 2, exhibit quite a linear 
behavior for (romTr)^ > 2, i.e., TO,r > 500 MeV; no ap- 
preciable curvature characteristic of the chiral logarithm 
is observed. 

One may suspect that pions in the simulation are too 
heavy to validate the use of ChPT. Another possibility 
may be the effect of explicit chiral symmetry breaking of 
the Wilson quark action at finite lattice spacings, a,s was 
discussed recently in the context of ChPT by [la| (see 
also llilll). Here we explore the more naturally look- 
ing, former possibility that the pion loop is suppressed 
for heavy pions and that the chiral logarithm manifests 
itself only for sufficiently small sea quark masses. The 
authors of proposed a model that incorporate such 
a behavior by introducing a hard cutoff regularization of 
the one-loop ChPT calculation. This amounts to the re- 
placement , m J In ml //^^ — > In m^ / {m^ where p. 
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FIG. 2: Chiral extrapolation of divided by the renormal- 
ization factor Za- The fits with the hard cutoff chiral loga- 
rithm are shown for = 300 (thin dashed curve), 500 (thick 
curve) and oo (dashed curve) MeV. 



is the scale of the hard cutoff, beyond which pion loop ef- 
fects are suppressed. This function has to be understood 
as a model when used above the cutoff \x. We may use 
this model to explore the possible range of uncertainties 
consistent with the lack of curvature in our data. 

Curves in Figure |21 illustrate the chiral extrapolation 
using the cutoff-logarithm plus a quadratic term. All 
curves are consistent with the lattice data, and they de- 
viate from the polynomial only in the small mass region. 
The fi = CO limit corresponds to the usual chiral log- 
arithm function plus a quadratic term, for which the 
curvature cancels between the logarithmic and quadratic 
terms in the region of simulations while giving a large ef- 
fect below (romTr)^ < 1. The other limit fi — MeV cor- 
responding to the polynomial fit. The variation with the 
parameter fi is taken as uncertainties in the chiral extrap- 
olation within this model. This shows that the value ob- 
tained with the polynomial fit (// = MeV), 147(3) MeV 
(here the errors are statistical only), may be affected by 
a 'hidden' chiral logarithm, leading to 128(2) MeV if the 
effect is maximal (/i = oo MeV). 

A similar analysis can be made for the heavy-light de- 
cay constant we have discussed above. The fits are shown 
in Figure^for fi — 300 and 500 MeV (thin dotted curves) 
as well as for oo MeV (dashed curve). The effect of the 
chiral logarithm can be —11% for /s, if we take fi — oo 
as an extreme case. While this case is probably unre- 
alistic, as it implies the validity of ChPT at very large 
mass scales, we take this as a conservative estimate of 
the systematic error, giving the lower limit for fp. Other 
functional forms may also be adopted (see, e.g., |20j), but 
such models are expected to give numerically similar re- 
sults in so far as the model is constrained by lattice data 
in the heavy pion mass region and by ChPT in the light 
pion mass region. 



The effect of the chiral logarithm is small for Jb^ , since 
the particle circulating the loop is kaon or eta. The for- 
mula in the partially quenched QCD is given in (ij . The 
chiral extrapolation is shown in Figure ^ with the lines 
for two extreme cases ii — and oo MeV. The difference 
between the two is only 1%. 

To quote our results we take the central value from the 
polynomial fit and include the variation in the presence 
of the chiral logarithm as an error. We obtain 



191(10)(ti"9)(12)(-) MeV, 
215(9)(t°)(13)(i^) MeV, 



1.13(3)(±io2)(2)(i3 



(2) 
(3) 

(4) 



where the first error is statistical, the second is the un- 
certainty from the chiral extrapolation and the other two 
are systematic errors explained in what follows. The 
error from the chiral extrapolation is one sided, since 
the polynomial fit is taken as our central value. The 
systematic error given in the third parenthesis is those 
arising from the finite lattice spacing (truncation of the 
actions and currents, and their perturbative matching) 
and the truncation of terms higher order in 1/ mg in the 
NRQCD action. An order estimate of the truncation 
errors shows that the most important contributions are 
0(A2jCD/"ifc) 4% and 0{al) ~ 4%. We add these er- 
rors by quadratures together with other possible errors. 
In , it is shown for the quenched lattice that such esti- 
mates correctly describe the error of finite lattice spacing 
(see also fo^' Bb). The errors we obtained is consis- 
tent with those in the quenched case at the comparable 
lattice spacing. The errors in the last parenthesis repre- 
sent the ambiguity in the determination of the strange 
quark mass, for which we adopt the value with the K 
mass as the central value. 

For the B parameter, ChPT predicts -(1 - 3g^)/2 
for the coefficient of the chiral log term instead of 
3(l + 3g2)/4 in ^ Therefore, the effect of the chiral 
logarithm is negligible in practice. For Bb^ there is no 
chiral logarithm as a function of sea quark mass. 

Figure 13 shows the chiral extrapolation of Bb (fJ^b) at 
Mb = iT^b (= 4.8 GeV) and the fits without the chiral 
logarithm. The triangles show the quenched results. The 
sea quark effect is small for this quantity. 

Our unquenched results obtained with a linear chiral 
extrapolation are 



BBjmb) 

BBA^b) 

Bb^ 
Br^ 



0.836(27)(i27)(56)(- 
0.850(22)(ii«)(57)(tg), 

1.017(16)(t5o')(17)(l^). 



(5) 
(6) 

(7) 



The meaning of errors is the same as for fs, except for 
the second one, i.e., those associated with the chiral ex- 
trapolation: we take the central value from the linear fit 
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FIG. 3: Chiral extrapolation of Bb^ (filled circle) and Bbs 
(open square). The four-quark operator is perturbatively 
renormalized, and the "method I" is employed for the trun- 
cation of higher order corrections (in as and l/mt) as an 
example [TH|. 



and put the difference from the linear plus quadratic fit 
as the systematic error. 

The amplitude of neutral B meson mixing is pro- 
portional to fgBg. Using the conventionally adopted 
renormalization-scale independent definition Bg, which 

(5) _ 



is related to BB{mb) as Bb 
225 MeV, we find 



1.528Bs(mf,) for A: 



MS 



Y Bb^ 



215(ll)(t23)(15)(-) MeV, (8) 
245(10)(t^)(17)(t^) MeV, (9) 



and for the SU(3) breaking ratio ^, 



(10) 



The chiral extrapolation gives the largest entry of sys- 
tematic errors for ^, as also suggested in [i^. Compared 
to the commonly assumed number in the phenomenolog- 

ical analysis, /s^-/^ = 230(28)(28) MeV {e.g., [H), 
where the second error is the quenching uncertainty, our 
central value of ||SJ) is slightly lower and the quenching 
error is eliminated. 

In conclusion we have obtained unquenched lattice es- 
timate of the B° - B° mixing matrix elements, including 
the decay constant, in a consistent set of simulations. Al- 
though the simulation is made at a relatively large mass 
of dynamical quarks, we explored the range of errors as- 
sociated with the chiral extrapolation: we expect that the 
true values of the matrix elements are within the range 
of indicated errors, even if the chiral logarithm would 
become manifest at a small quark mass. 
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